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abstract 

A computer program has been developed which generates Feynman 
graphs automatically for scattering and decay processes in non-Abelian 
gauge theory of high-energy physics. A new acceleration method is 
presented for both generating and eliminating graphs. This method 
has been shown to work quite efficiently for any order of coupling 
constants in any kind of theoretical model. A utility program is also 
available for drawing generated graphs. These programs consist of 
the most basic parts of the GRACE system, which is now used to 
automatically calculate tree and one- loop processes. 
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1 Introduction 



In regard to recent projects concerning experimental high-energy physics with 
increasing available energy, accurate theoretical analysis has been required 
for increasing the number of scattering processes. Since electro-weak theory 
combined with QCD, known as the standard model, is considered to be a ba- 
sic theory among various theoretical models, exact perturbative calculations 
within the framework of this model are considered to be the standard for the- 
oretical predictions. Because of the complexity of interactions, the number 
of Feynman graphs for one process becomes much greater than that in QED. 
There appear several tens to hundreds of graphs in the tree process and sev- 
eral hundreds to thousands in the one-loop process. In addition, in contrast 
to pure QCD, one cannot ignore various mass parameters which play impor- 
tant roles. These situations make a theoretical calculation of the amplitude 
more difficult and tedious. At present, it seems that the amount of labor 
necessary for exact calculations has almost reached the limit of hand calcula- 
tions. It is thus natural to utilize computers to carry out such work. Several 
groups have started independently to develop computer systems which au- 
tomate the perturbative calculations in the standard model 0, |2|, [| f|, [J. 
Automatic calculations of the tree processes and some part of the one-loop 
processes have already been achieved by these systems. 

An automatic calculating system starts with the generation of Feynman 
graphs for a given physical process. Feynman-graph generation for all orders 
in QED is not difficult, since the problem can be reduced to a problem equiv- 
alent to tree-graph generation^, [/J. For electro-weak theory, or even in the 
(j) 3 model, however, the situation is different. They include self- interactions 
of particles, whose vertices are symmetric under exchanging interacting par- 
ticles. This symmetry of the vertices results in a complicated structure of 
the internal symmetry in a Feynman graph. 

It is not difficult to generate a sufficient set of Feynman graphs by sim- 
ply connecting the vertices. The problem, however, is to eliminate duplicated 
graphs. It is hard to analyze the structure of the internal symmetry of graphs 
so as to avoid duplicated graphs controlling the graph generation process. 
Usually, a newly generated graph is discarded when it is found to be topo- 
logically equivalent to one already generated. Although explicit comparisons 
of graphs are, of course, possible, any known algorithm consumes time, which 
increases as an exponential of the number of vertices. A graph comparison 
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is generally not considered in the category of polynomial time complexity^ . 
There is another problem in a explicit comparison method. It requires to save 
all generated graphs in order to recall graphs for comparisons. In a practical 
problem, the list of generated graphs grows too long to keep in the main 
memory of a computer. So the list is recorded on a secondary memory, usu- 
ally on disk space. Frequent access to a secondary memory makes program 
very slow. This problem is solved by an orderly algorithm^. An orderly 
algorithm judges whether a newly generated graph is necessary to be kept or 
not, looking only at the graph without recalling already generated graphs. 
Based on this algorithm, a Feynman-graph generator for a wide variety of 
models was developed by P.Nogueira||TO| . 



The aim of this paper is to present a computer program for Feynman- 
graph generation that is sufficiently fast for practical use. For this purpose, 
we have developed a new method of graph generation which accelerates the 
orderly algorithm. We classify vertices in such a way that there is no topo- 
logically equivalent vertices in different classes. With this classification, we 
are able to reduce the number of graphs to which the orderly algorithm is 
applied. The classification method of vertices which we use is an empirical 
one developed by graph theorists in order to determine the equivalence of 
two given graphs ]ll|. We have modified this method so as to be suitable for 
graph generation. 

In order to make the problem simpler, we assume the following conditions: 

1. Vacuum-to- vacuum graphs are not considered. 

2. External particles of the graphs are assumed to be topologically differ- 
ent. 

The first condition stems from the fact that the purpose of our automatic 
system is to calculate physical processes, namely the scattering and decay of 
particles. With the second condition, we have a one-to-one correspondence 
between the Feynman graphs and the Feynman amplitudes. One can gener- 
alize the method to loosen these restrictions without any essential difficulty, 
since the orderly and classification algorithms are still applicable. However, 
it must make the program more complicated. This generalization is discussed 
in section [6|. 

In the next section we introduce our basic methods. An explanation 
concerning the orderly algorithm is also given. 
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In order to accelerate orderly algorithm, we classify vertices in section 
H A pre-selection rule of generated graph, which is necessary to make the 
orderly algorithm consistent with vertex classification, is presented in this 
section. 

The method of graph generation is described in section [|. The vertex- 
classification algorithm has been combined in order to eliminate graph in 
intermediate steps. 

In section |5] we introduce our implementation of algorithms and utility 
programs. 

A summary and comments are given about the physical processes calcu- 
lated by the GRACE system up to now. A generalization of our program 
without the two conditions mentioned above is also discussed. 

Appendices [A|, ||, |C| and |D| comprise a brief manual of the programs and 
file formats. 

2 Basic method 

A graph consists of a finite number of nodes and edges. A node is either an 
external particle or a vertex. An edge is either a propagator or a connection 
between an external particle and a vertex. Throughout this paper, N denotes 
a set of nodes and E C {(u,v)\u,v G N} does a set of edges. A graph 
G = (N, E) is defined mathematically as a pair of these sets. The degree 
deg(v) of node v is defined as the number of edges attached to the node v: 

deg(v) = \{(u,v)eE}\ (1) 

The degree of an external particle is 1. An edge has several attributes in a 
Feynman graph, such as the number of multiple connections between nodes 
and the name of a particle of a propagator. In order to make our discussion 
simpler, we drop these attributes in the following description. A node is 
labeled by a positive integer, with which we identify a node. We regard a set 
of integers {1, \N\} as the set of nodes N. 

Let us consider a graph G = (N, E) and a permutation p mapping the 
set of nodes N onto itself. The permutation p is generalized to act on the 
set of edges by defining p(u,v) = (jm,pv) for (u,v) G E, u,v G N, and is 
generalized to act on the graph G as pG = p(N,E) = (pN.pE) = (N,pE). 
Two graphs, G = (N, E) and G' = (N, E') , are isomorphic when there exists 
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such a permutation p that maps the graph G to G' = pG. This condition is 
equivalent to 



Isomorphic graphs are different each other only in the way of labeling nodes. 

A permutation p is an automorphism of a graph G when pG = G. The set 
of all automorphisms of a graph G forms a group Tq called an automorphism 
group. 

A graph can be expressed as a matrix, called adjacency matrix, whose 
element carries the number of connections between nodes % and j. We 
show an example of isomorphic graphs and their adjacency matrix in Fig. |T[ 
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Figure 1: Example of adjacency matrix. 
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In this and the next sections we consider a simple method of graph gen- 
eration. First, necessary number of nodes are generated and then the nodes 
are connected in all possible ways. This method creates all possible graphs. 
Topologically equivalent graphs may appear several times. In order to elimi- 
nate duplicate graphs, we use a systematic method, an orderly algorithm^, 
which judges a newly generated graph being necessary or not without recall- 
ing any of already generated graphs. 

An orderly algorithm compares graphs through values of a function, called 
coding, which maps a graph to an integer in a way that different isomorphic 
graphs are mapped to different values. Such a function is easily realized, for 
example, in regarding the elements of the adjacency matrix as a sequence of 
digits of an integer 0. 

In order to pick up a representative from a set of isomorphic graphs, An 
orderly algorithm selects a graph G when it satisfies the following condition 
with coding /: 

f(G) = max f{pG), (3) 

pe5|jv| 

where Sm\ is the symmetry group of the set of nodes N (the set of all per- 
mutations acting on N). A graph G is indirectly compared with graph pG 
through their values of the coding /. 

Combining this selection condition with graph generation process, an or- 
derly algorithm generates graphs in the following procedure: 

1 Generate all the possible graphs. 

2 Apply all possible permutations to each graph. 

3 Discard a graph when it has a smaller coding value than ones of per- 
muted graphs. 

This method requires to keep only two graphs, a tested graph and its 
permuted one, which are small enough to keep in the main memory of a 
computer. 

When a permutation p satisfies selection condition (Q), p is an element of 
automorphism group Tq of G. One can construct Tq explicitly in collecting 
all of such permutations. The symmetric factor |r<y| necessary for calculating 
Feynman amplitudes is evaluated in this way. 
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A simple application of this algorithm requires comparisons among graphs 
corresponding to 0(|iV|!) permutations for each of 0(\ N\\) isomorphic graphs. 
In order to decrease the number of comparisons in the orderly algorithm, we 
consider to replace S\m in (|3[) by its subgroup. 



3 Vertex classification 

We consider a classification {iVj} (UjiVj = N, Nj fl Nj = for i ^ j) of the 
set of nodes N such that any two nodes in different classes are topologically 
different. We do not require that nodes in a class are topologically equivalent. 
In other words, each class of nodes is a direct sum of the orbits in (an 
orbit of a node v in Tc is the set {it|3p G Tc,pv = it}). We call such a 
classification a consistent classification. With this classification of nodes, we 
construct a group 

S({N t }) = S\ Nl \®S m ®... (4) 

where is the symmetry group acting on the set iVj. We replace the group 
S\n\ in condition (H) by this subgroup S({Ni}). This replacement reduces 
the number of permutations from |iV|! to |iVx|! x |AT 2 |! x .... Corresponding 
to this limitation of comparisons in the orderly algorithm, one must limit 
the set of generated graphs. We present the way of limiting graphs after 
describing our method of node classification. 

We adopt a method of node classification used in graph theory [p]. First, 



nodes are classified into {N- ^} by simple topological properties. Since we 
consider that external particles are topologically different from each other, 
they are put into different classes. The i-th external particle is put to the 
i-th class. Other nodes are classified by their degree such that: 

u, v G <^=^ deg(it) = deg(t>) 
u G , v G Nf\deg(u) < deg(v) ^i<j (5) 

We call this initial classification {N^} a primitive classification. Primitive 
classifications {N- ^} and {M- ^} of isomorphic graphs G and pG, respec- 
tively, satisfy: 

pN^ = Mf } for all i, (6) 
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since relabeling of nodes does not change topological properties of nodes such 
as deg(pv) = deg(f). 

Starting from the primitive classification, the classes {N-^} are refined 
to {N^ k+1 ^} iteratively. We consider a vector a v [i},i = 0,1,... for a node 



v G Nj k \ The zeroth element o„[0] of the vector keeps its current class 
number, and the z-th element keeps the number of edges which connect v 
with nodes in the i-th class N^: 



a v [Q] = j, for v E N, 



(k) 



a v [i\ = \{(u,v) G E\ue N^ h) }\ forz>0. (7) 

The nodes are classified by the values of these vectors. New classes are num- 
bered in increasing order, which is evaluated in the lexicographical ordering 
"-<" of vectors. This process is repeated until classes can no more be re- 
fined. Thus the nodes are classified only by their topological properties. We 
summarize this refining method of classification as follows: 

refine ({7Vf } }) 
{ 

for all (class j) { 

for all (v G A/f } ) { 
a v [0] = j; 
for all (class i) 

a v \i] = \{(v,w)eE\weNt ) }\; 

} 

} 

(construct {N- k+1 ^} such that 

a v = a w ,ve N^ k+1 \w G +1) i = j ; 
a v ^a w ,ve N^ k+1 \w G A/f +1) =>• % < j ; 

); 

return {iVf +1) }; 

} 

We call a sequence of this classification of nodes {N-^} ( k — 0, ...,K ) a 
refinement sequence. 

We show the following properties for a refinement sequence. 
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Property 3.1 For any N^ k+1 ^ there exists an input class Nj® includes N^ 1 ' : 

N (k+i) c N (k)_ (g) 

This property states that N- k+1 ^ is a refined class of Nj. 

Property 3.2 Ordering of classes are kept by the refinement procedure: 

Vi Vj [ivf +1) C , Nf +1) c JVg) , • < j =► I < m) . (9) 

This property is evident from the fact that the algorithm renumbers the new 
classes according to the lexicographical ordering of the vector a v , which keeps 
the old class number at the first element. 
Thus classes are divided to: 

N (k) = iV f+i) uA r(^) u ... u ^+i) 

N W = N^UN^U...UN^ (10) 



Property 3.3 When nodes of two isomorphic graphs G = (N, E G ) and H = 

pG = (N,E H ), p G S\m, are classified into {N^} and {M^}, respectively, 
the following relation holds for all k: 



Vi \pN$ k) = M\ k) ] Vj [ P N^ k+1) = Mf +1) ] (11) 

Proof. Let us take a node v G N- k \ We consider vectors a° and a^ v in 
eq. (0) which are used to construct {Nj k } and {Mj fc+1 ^} from {iV"j } and 
{Mj k ^}, respectively Their zero-th elements are equal a G [0] = a^[0] = i 
since pv G . We obtain for I > 0: 

a G v [l] = \{(v,w)eE G \weNl k) }\ 

— \{(P V ,P W ) G pE G \pw G pNi k ^}\ 

= \{(pv,pw) G E H \pw G Mf } }| (12) 

= <M 
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Thus equation a„ = a^ v is hold for any node v. Since the numbering of 
refined classes are determined only by these vectors, it is easy to see 

v e ivf +1) <=> P v e Mf +l) , (13) 

which implies pNj k+1 ^ = Mj k+1 \ □. 

Since pN- ^ = M- ^ is satisfied for all i in eq. @, the above property 
result in: 

pN$ k) = M\ k) for all i and k. (14) 

Especially, the length of refining sequences of two isomorphic graphs are 
equal. 

Now we consider the way of limiting graphs to be applied to the orderly 
algorithm. We impose the following pre-selection condition for the classifica- 
tion {iVj} of nodes of a generated graph: 

Vi Vj Vu e N t Vv e Nj [i < j => u < v], (15) 

where the ordering of the nodes is evaluated in terms of the labeling numbers 
of the nodes. Classes satisfying this condition are expressed as: 

n, = {l,...,!^!} 

N 2 = {\N 1 \ + 1,...,\N 1 \ + \N 2 \} (16) 



We consider this condition not only for the final classification but also clas- 
sifications appearing in a refinement sequence. 

Property 3.4 Classification {N^} is necessary to satisfy condition ftTEj) 
for {N- h+1 ^} to satisfy the same condition. 

Proof. Let {N^ k+1 ^} satisfy the condition ([15]). A class Nj k+1 ^ is expressed 
as eq. ([16]): 

Nf +1) = {sf +1) + 1, sf +1) + 2, sf + \% (17) 
where sf +1) = Ekj Wf +1) \. 
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Let a class is decomposed to b{ classes of {N- K+L> }. From eq. (|10D, 



(fc+i)- 



we obtain, 



N. 



(k) -= U U ... U Ni k+1 \ 



where Cj = J2j<ibj- The following relation holds: 



Z<Cj + l 

i=Ci_l + l 



(*+l)l 



(fc+1) , I AJW 



(fc) 



ffc+1) 



Combining with c\ — and Sci+i = = 0, we obtain 



+ s> - s) 

.(fc) 



Equations flU), (0) and © lead 



iV. 



(fc) 



l ,6, r Q (fc+!) , i «.(*+!) , o ,(*+!) 1 
f„(*+l) , i ,(*+!) 1 



This expression is in the form of eq. (|T6|). Thus {iV^ } satisfies condition 

© □. 

With this property, we can eliminate a graph before reaching to the end 
of the refinement sequence, when an intermediate classification in a classifi- 
cation sequence is found not to satisfy condition (|T5|). So we can decrease 
the number of refinement steps for graphs which are to be discarded. 

We show that graph selection using this condition (|T5|) leaves enough 
graphs. 

Property 3.5 For each graph G, there exists such a graph that is automor- 
phic to G and satisfies condition fill ). 



(18) 



(19) 



(20) 



(21) 
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Proof. We consider the case when {A^} does not satisfy condition (jl5|). We 
define a permutation p of nodes by the following procedure: 

void renumber () 
{ 

count = 1 ; 

for (i = 1; i <= (The number of classes) ; i++){ 
for all (v e { 
(let pv = count) ; 
count = count + 1 ; 

} 

} 

} 

It is easy to see that classification {pNi} is expressed as eq. ([16]). So the 
graph pG satisfies condition (|T5|) □. 

We select graphs by the following condition instead of eq. (Q): 

f(G) = max f(pG) (22) 

Our modified orderly algorithm becomes: 

1 Generate all the possible graphs. 

2 Select graphs by condition ([i~5]). 

3 Apply all permutations in S({Ni}) to each graph. 

4 Discard a graph when it has a smaller coding value than ones of per- 
muted graphs. 



In order to show that two selection rules flTq ) and (|22|) are consistent each 
other, we prove the following relation for a graph G satisfying ([15]): 

{pG\p E S ]N{ ,pG satisfies condition (15)} = {pG\p E S({Ni})}. (23) 

Property 3.6 When nodes of a graphs G are classified into {A^} satisfying 
condition ftT3{ ) and a graph H = pG is isomorphic to G with permutation p, 
necessary and sufficient condition that a graph H satisfies condition j^T^ ) is 
that the permutation p is an element of the subgroup S({Ni}). 
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Proof. Let the nodes of isomorphic graph H = pG be classified into 
{Mj}. From eq. (0), a class iVj of G is mapped by p to the class Mi = pNi 
with the same index. 

When a permutation p is an element of 5({A^}), each class iVj is kept 
invariant under p, that is, N{ = pN^ = Mi. Since {Aj} satisfies condition 
(p~5|) , {Mj} satisfies the same condition. 

Inversely, we assume that the graph H satisfies the condition flilf) . Since 
|iVi| = |piV"i| = | M»| holds for all i, it is easy to see from expression (|I6|) that 
Ni = Mi = pNi is satisfied for all z. This implies that the permutation p 
keeps {Aj} invariant and that p is an element of S({Ni}) □. 

We have first confirmed that necessary graphs are kept in step 2 of our 
modified algorithm. Then surviving isomorphic graphs are transformed each 
other by an element of subgroup ^({Aj}). They are compared through the 
values of coding in step 4 and duplicated graphs are properly eliminated. 

In the next section we combine step 1 and 2 in order to decrease the 
number of generated graphs. 

4 Graph generation and selection 

Here we describe how to generate graphs in combination with elimination of 
graphs. 

The first step of graph generation is to prepare nodes. A node is to be 
connected to several edges. We consider such an imaginary object, let us call 
it a leg, that is a place in a node where an edge is to be tied to. A fixed 
number of legs is assigned to each node at the beginning. The number of legs 
becomes equal to the degree of the node when connection process terminates. 

The next step is to connect nodes by edges. The nodes are connected 
iteratively starting from one fixed external particle, which we call the root 
of the graph. In order to make our way of connections systematic, we define 
level of node v as the distance of the node from the root measured by the 
minimum number of edges among paths connecting v and the root. In this 
step there appear intermediate configurations of graphs, in which only a part 
of legs of nodes are tied with edges and some nodes may remain isolated from 
other nodes. We consider the level of an isolated node is infinity. 

Nodes are connected in the following way. The root, with one leg, is first 
connected to another node. The connected node becomes the only one node 
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at the first level. Then this node is connected to other nodes until all of its 
legs are connected. These connections define the set of nodes in the second 
level. Then, the all legs of all nodes in the second level are connected to 
others. In this way, nodes are connected in increasing order of the value of 
level. Once a level of a node has a finite value, its value is not changed. This 
process proceeds in a recursive way so as to exhaust all possibilities of the 
connections. A skeleton description of the algorithms is given as follows: 

void gsconnO 
{ 

node ns , nt ; 



ns = (find node at the lowest level with free legs) ; 
if (no more such nodes) { 

if (the graph is a connected graph) { 

if (the graph is accepted by the orderly algorithm) 
(a new graph is obtained) ; 

} 

} else { 

for all (node nt with free legs) { 
(connect ns to nt) ; 
gsconnO ; 

(disconnect ns from nt) ; 

} 

} 

} 

For efficient graph generation, it is important to eliminate unnecessary 
intermediate graphs in as earlier stage as possible, since a large number of 
graphs may be produced in exhausting all possible additional connection to 



the intermediate graph. We eliminate graphs with condition (15) not only in 
the final form of the generated graphs, but also in intermediate graphs. 

One must notice that the classification of an intermediate graphs is not 
always consistent with that of the final form of the graphs. For example, 
we consider a final form of the graph of Fig. la and its intermediate graph 
Fig. lb. The set of classes of nodes of Fig. la is 

{(1), (2), (3), (4, 5), (6)}, 
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Figure 2: Example of a two- loop graph and its intermediate configurations. 

while that of Fig. lb is 

{(1), (2), (3), (4), (5), (6)}. 

This example shows that nodes 4 and 5 are topologically equivalent in the 
final form of the graph, although they are not always equivalent along the 
way of graph generation. When one applies condition fll5|) to such graphs, 
one looses necessary graphs. It is necessary to apply the condition flT5D only 
to the consistent classification. 

We can recover a consistent classification by considering an intermediate 
graph = (N, E^) at the time when all of the legs of all of the nodes in 
level I are just connected. These intermediate graphs form a finite series, 

q[o] c G [i] c _ c Q[m-1] = Gj 

where m is the maximum level of the nodes in graph G. Fig.lc represents an 
example of G^. The set of edges is expressed as: 

E® = {(u, v) E E\ min(level(u), level(w)) < I}. (24) 

We show the following property for the automorphism group r Gra of CW: 
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Property 4.1 The following relation holds for all I: 

r G [i+i] C T G [i]. (25) 

Proof. If this is not satisfied, there exists a permutation p e r^in-i] which 
is not an element of T G [i]. This means that there exists an edge (it, v) G 

E [i] 

and (pu,pv) ^ E^. On the other hand, since p is an automorphism in T G \i+\\ 
and level (pu) = level (it), we get: 

min(level(pu), \evel(pv )) = min (level (it), level (v)) < I, (26) 

which implies (pu,pv) E E^ □. 

Relation (|25|) implies Tg = r^m-i] C r G [i]. Thus an orbit of r G [i] is a 
direct sum of the orbits of T^. Since a class of nodes constructed for is a 
direct sum of orbits of the automorphism group T G [t\ , the class is a direct sum 
of the orbits of T G . The classification of nodes in is therefore consistent 
with the classification in G. 

We can now eliminate irrelevant intermediate graphs with condition (|15|) 
applying to G^. Furthermore, it is possible to use the classification of nodes 
in 

GrP-i] 

as an input for constructing the classification in in our class 
refinement algorithm. This method decrease the total number of class 
refinement. The final form of our algorithms is as follows: 

void gsconn(int lvl) 

{ 

node ns , nt ; 

if (all legs of the nodes in the level lvl have been 
connected) { 

(refine classes of the nodes) ; 

if ( ! (class ordering condition fllSD is satisfied) ) 

return; 
lvl++; 

} 

ns = (find node at the level lvl with free legs) ; 
if (no more connectable legs) { 

if (the graph is connected one) { 

if (the graph is accepted by the orderly algorithm) 
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(a new graph is obtained) ; 

} 

} else { 

for all (node nt with free legs) { 
(connect ns to nt) ; 
gsconnO ; 

(disconnect ns from nt) ; 

} 

} 

} 



5 Implementation 

We have implemented the method described above as a computer program 
written in C language. The program generates Feynman graphs in which 
particles are assigned to propagators. The particles of the propagators are 
determined in accordance with a table of particles and interactions defined by 
the users. It is possible that duplicated graphs are produced in this particle 
assignment process. They are eliminated again by the orderly algorithm with 
extended coding of graphs with particle attributes on the edges. 
The program has the following options for graph generation: 

1. To pick up only one-particle irreducible graphs ]13 . 

2. Not to generate graphs with a self-energy part at an external particle 
line. 

3. Not to generate graphs with a self-energy part at an internal particle 
line. 

4. To generate graphs with counter terms for renormalizing the theory. 
Counter terms are automatically generated when their interactions are 
of the same form as tree vertices. Other counter terms can be added 
by users (see appendix y). 

5. To generate a skeleton graphs, in which looped one-particle irreducible 
subgraphs are considered as to be blobs. The generated graphs are of 
the tree type, but with blob vertices. 
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6. Not to assign particles to the internal lines (only topology). 



The way to specify these options is described in appendices [X| and |B|. 

We also provide a graph-drawing facility on the X-window system and on 
PostScript files. We show an example of output figures in Fig.0. 
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Figure 3: Example of drawn graphs. 



The total number of graphs depends on the conditions of the graph se- 
lection. In special cases, they are analytically calculable. In the 4 model 
including <j) 3 interaction, the method of enumerating connected graphs is de- 
veloped in a graph-theoretical method W\ [L5[ . When there appear different 



kinds of particles in the model, a weighted sum by the symmetric factor 
can be calculated in zero- dimensional field theory ||16|| . These numbers of 
graphs are calculated and listed in ref.[|l(|. We have checked our program 
by comparing with these numbers. Since for up to 5-loop graphs the orderly 
algorithm is not necessarily required to generate unique graphs fl2| , we have 
also checked the number of graphs for 6- loop tadpole, which includes 90156 
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graphs[]. The number of one-particle irreducible graphs is checked using a 
recursion formula for one- and two-loop graphs. 
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11460 


21.3 
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2485 


121520 


41.1 






2485 


3.3 



Table 1: Efficiency of acceleration in the 4 model including 3 interaction. 

The performance of the program is shown in table |I| in the case of con- 
nected graphs in the same model. Columns A, B and C represent the num- 
ber of external particles, the number of loops and the number of graphs, 
respectively. The theoretical prediction of the number of graphs is calcu- 
lated by using a program written in formula manipulate language REDUCE, 
which takes more cpu-time than our graph-generation program. Columns 

3 The corresponding value in ]l2t is not correct. 
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D, E and F correspond to simple application of the orderly algorithm (D), 
application of condition without (E) and with (F) elimination of inter- 
mediate graphs. Sub-columns G and H are the number of graphs applied to 
the orderly algorithm eq. (^) or (|22]) and cpu-time in seconds measured on 
SUN-IPX with "-a -o /dev/null" options (described in appendix |A|). This 
result shows that elimination of intermediate graphs is effective enough; it 
reduces cpu-time in several order of magnitudes. 



6 Summary and comments 

We have developed a method of Feynman graph generation, which accelerates 
the orderly algorithm, based on a node classification method and systematic 
elimination of intermediate graphs. This method is implemented as a com- 
puter program. The order of the coupling constants of a physical process 
is not limited. Graphs can be generated based on a user-defined model; in 
this paper we adopt electro-weak theory combined with QCD as an example. 
The program has several options of graph selection, including to generate 
renormalization counter terms. Although the asymptotic behavior of the 
execution time is proportional to the factorial of the number of nodes, the 
program is sufficiently fast for practical use. 

This program is a part of the GRACE system, which automates com- 
putation of tree and one-loop scattering processes. With this system, the 
cross sections including one- loop corrections in the processes e + e~ — > HZ 
and e + e _ — > ti, and one-loop gluon corrections to the process e + e~ — > qqj 
are automatically calculated]^]. For the tree case, the cross sections of the 
scattering processes with up to five final particles are calculated as shown in 
table 0. The number of graphs is counted in covariant and Feynman gauges 
for tree and one-loop processes, respectively. In both cases, the interactions 
among a Higgs particle and two light fermions are neglected. Integration 



over phase space is calculated using the integration package BASES fl7|l. 
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Process 


#Graph 


Reference 
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l|a] 


-» v e u e Z° 


9 


l 


l|b] 


— > v e 9 e H 


2 




18|b] 


-> e+e-H 


2 




lgjb] 




4 


1 


lsjb] 


-> VF+l^-ff 


11 




18b] 


-> z° z° z° 


9 


1 


l|b] 




20 


l 


lSjb] 


-»ttz° 


9 




l§b] 




6 




18b] 


—> H H Z° 


6 


l 


l|b] 


-> 777 


9* 




77 — > e+e~Z 


6 




7 e -> eW+W~ 


18 




99 -> qqj 


18 





Table 2a Tree processes with the final 3-body 



Process 


#Graph 


Reference 


e+e" -> v e v e W+W~ 


60 




18b] 




36 






-> e+e-W + W~ 


114 




18b] 


-» v e v e Z a Z a 


57 




lip] 


-> e+v e W-Z° 


88 




18b] 


-^e+e-Z°Z° 


86 




18b] 


— > e+e-77 


80 




lid] 




40 




l|d] 


— > v e D e bb 


21 




18b] 




11 




18b] 


— » e + v e tb 


21 




is|f] 


_> V^+VF-77 


138 




— » erv e ud 


24 


1 


lis] 


— > udud 


69 




-> gg#7 


21 




— » e + e~bb 


50 




7e -> v e W-HH 


40 





Table 2b Tree processes with the final 4-body 
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Process 


#Graph 


Reference 




71 




— > e~v e udr( 


142 





Table 2c Tree processes with the final 5-body 



Process 


# Graph 


Reference 




89 






-> tt 


50 


1 


1 



Table 2d One-loop processes with the final 2-body 



Process 


#Graph 


Reference 


e + e~ — > gg7 


12 





Ghion correction with only final radiations. 



Table 2e One-loop processes with the final 3-body 



Table 2: Scattering processes calculated by the GRACE system. 
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In the following we discuss how to loosen the restrictions for graph gener- 
ation assumed in the introduction. The second restriction regarding identical 
external particles can be removed by changing the primitive classification so 
as to put identical particles in the same class. At the same time, the root 
node should be replaced by a class of nodes. 

The first restriction of excluding vacuum-to- vacuum graphs is complicated 
to loosen. It is possible that the primitive classification contains only one 
class, such as in vacuum-to- vacuum graphs in the <f> 3 model. In this case, node 
classification and acceleration methods do not work anymore. However, one 
can artificially fix one node as the root of the graph. 

We consider a set of graphs A generated by our method with fixing node 
1 among n = \N\ nodes and selected by condition flip]). Let us also consider 
the true set of vacuum-to-vacuum graphs T. For any graph GeT, one can 
construct a set of isomorphic graphs Iq = {pG\p G S n _i}, where 5 n _i is 
the symmetric group acting on the set of the nodes {2, • • • , n}. We consider 
to select a representative graph from Iq by applying condition (|15|) after 
classifying the nodes with fixing node 1. The selected graph belongs also to 
A, since A is constructed by the same selection rule applied to all possible 
graphs. The set of graphs A thus contains at least an isomorphic graph to 
any graph of T, though some graphs in A are duplicated. 

The duplication of graphs is to be tested by the original orderly algorithm 
(H). Since we have already eliminated some of the duplicated graphs by 
condition (|15|), we must limit the set of graphs {pG \ G G A, p G P} by 
the same condition before comparing the values of the coding. The orderly 
algorithm is then changed to the following form: 

compare (GO 

{ accept = True; 

for all (pG, pE S\n\) { 
H = pG; 

(classify node of H with fixing node 1) ; 
if (condition (|15|) is satisfied) { 

if (coding (G) < coding (H)) { 

accept = False; 

break; 

} 

} 
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if (accept) 

(accept graph G) ; 

} 
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A Installation and execution 



The source code is available by anonymous ftp from ftp.kek.jp in the 
directory kek/minami/grc. It has two programs: one generates Feynman 
graphs; the other draws generated graphs, which work on UNIX with the 
X- Window system. 

The graph-generation program grc reads two files, a process file and a 
model file. It then creates an output file named "out.grf", which keeps 
the information about the generated graphs. The format of the input files is 
described in appendices |B] and [Cj 

The graph-drawing program grcdraw reads the model and the output 
file of the generation program, and then draws the graphs on the display 
with the X- Window system. This program is not so intelligent as to dis- 
play graphic objects beautifully; however, it is still useful for checking the 
generated graphs. This program has several subcommands. 

The procedure of installation is: 

1. Uncompress and expand the source code file. 
In a new directory, run the following command: 

zcat grc. tar. Z I tar xvf - 

2. Editing Makefile. 

The file Makefile in the sub-directory src controls the way of com- 
piling source code. The path-names of include files and libraries of 
Xlib should be changed in accordance with the system configuration. 
If Athena widgets (included in the standard distribution of X- Window 
system) is installed in your system, use the line in the Makefile: 

XCOPT = -DTOOLKIT -DX11 

If not, use 

XCOPT = -DX11 

3. Compilation. 

In the sub-directory src, run make command. 
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4. Testing the Feynman-graph generator. 

In the sub-directory sample, run the following command 

. . /src/grc 

The program generates 4 tree graphs for the e + e~ — ► W + W~ process, 
316 one-loop graphs of the same process and 28 tree graphs for e + e~ — > 

5. Testing the Feynman-graph drawer. 

In the same sub-directory sample, run the following command 

. . /src/grcdraw -h 
. . /src/grcdraw 

Try to type "f", "b", "n", "p'\ "g" V- 

If you have installed with Athena widgets, buttons will appear on the 
display. 



Feynman-graph generator 

The command syntax of grc command is 

grc [options] [.process-file-name] 

The process-file-name is an input file specifying the physical process and 
options, the format of which is described in appendix |B|. The default process 
file name is "in. pre". 

The command line options of the grc command are: 

• -a 

Skip particle assignment. 

• -c 

Do not generate an output file, but only count the number of graphs. 

• -o output-file-name 

Specify the output file name. The default output file is "out.grf". 
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• -h 

Print help message. 
Also try with the help option: 

. . /src/grc -h 

Feynman-graph drawer 

The command syntax of the grcdraw command is 

grcdraw [options] [graph-file-name] 

The graph-file-name is an output file generated by the grc command, the 
format of which is described in appendix [D]. The default process file name is 
"out.grf". 

The command line options of the grcdraw command are: 

• -w number-of- graphs 

If this option is specified (number-of-graphs) 2 graphs appear on the 
display. 

• -h 

Print help message. 

Also try with the help option: 

. . /src/grcdraw -h 

This program has some subcommands for selecting graphs drawn on the 
display. They follows: 



q 

n 

P 
f 

b 

g 
s 

1 

t 

<number>j 
o 



quit. 

display next process 

display previous process 

display forward graphs 

display backward graphs 

scale up the size of graphs 

scale down the size of graphs 

display particle names of internal lines (on/off) 

display particle names and graph number (on/off) 

jump to the specified graph 

output the displayed graphs to PostScript file 
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When the program is installed with Athena widgets, such commands are 
also displayed as buttons. 

B Process file 

The graph-generation program works in accordance with the process file. 
The following example is a file specifying the one- loop and tree process of 
e + e" — > W + W~ and the tree process of e + e~ — > W + W~y. 

mmmmmxxxxmxmnmmm 

Model="all.mdl" ; 

Process ; 

ELWK={2, 4}; 

Initial={electron, positron}; 

Final ={W-plus, W-mimis}; 

Expand=Yes ; 

0PI=No; 
Pend; 
Process ; 

ELWK=3 ; 

Initial={electron, positron}; 
Final ={W-plus, W-minus .photon} ; 
Expand=Yes ; 
0PI=No; 
Pend; 

The lines beginning with "'/," are ignored as comment lines. 

The first non-comment line specifies the file which describes the model 
used for graph generation, which is described in appendix 0. 

Then follow descriptions of the processes. A block describing a process 
begins with the the line of "Process;" and ends with the line of "Pend;". 
In this block, descriptions are given for the order of the coupling constants, 
initial particles, final particles and some options. 

The name of the coupling constant ("ELWK" in the above example) is 
defined in the model file. When the value of the coupling constant is given as 
a list of numbers, the program generates graphs for each value of the coupling 
constants with the same external particles. 

The initial and final external particles are given as lists of particle names, 
which are defined in the model file. 
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The available options are : 

• OPI = Yes | No 

Generate one-particle irreducible graphs or not. 

• Expand = Yes I No 

Expand the looped part or generate looped parts as blobs. 

• Tadpole = Yes I No 

Generate tadpoles or not. 

• extself = Yes I No 

Generate the self-energy part at an external particle or not. 

• selfe = Yes I No 

Generate whether with the self-energy part or not at an internal particle 
line. 

• countert = Yes I No 

Generate renormalization counter terms or not. 
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C Model file 



Here, we describe the format of the model file. Our definitions of particles 
and vertices include not only the necessary information for graph generation, 
but also for amplitude generation. We show an example of this file: 

%%%%%%%%%%%%%%%%%%%7//.y.7.%%%%y.%%%%%%7.7.*/.7.%7. 

Order={ELWK, QCD}; 
%======================================= 

7. gauge bosons 

t 

Particle=W-plus ["W+"] ; Antiparticle=W-mimis ["W-"] ; 

PType=Vector ; Charge=l; Color=l; Mass=AMW; Width=AGW; 

MValue="80.22D0"; WValue="2 . 12D0" ; PCode=2; 
Pend; 



% 

% scalars 

% 

Particle=Higgs ["H"] ; Antiparticle=Particle ; 

PType=Scalar ; Charge=0; Color=l; Mass=AMH; Width=AGH; 

MValue="150.0D0"; WValue="0 . 0D0" ; PCode=3l; 
Pend; 

7. 



y> 

7. leptons 

% 

Particle=nu-e ["nue"] ; Antiparticle=nu-e-bar ["~nue"] ; 

PType=Fermion; Charge=0; Color=l; Mass=AMNE; Width=0; 

MValue="0.0D0" ; PCode=51; Massless; 
Pend; 



Particle=electron["e-"] ; Antiparticle=positron ["e+"] ; 

PType=Fermion; Charge=-1; Color=l; Mass=AMEL; Width=0; 

MValue="0.511D-3" ; PCode=55; 
Pend; 
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y 

7. quarks 

y 

Particle=u; Antiparticle=u-bar ["~u"] ; 

PType=Fermion; Charge=2/3; Color=3; Mass=AMUQ; Width=0; 

MValue="100.0D-3"; PCode=61; 
Pend; 



%======================================= 

7. vvv 

y 

Vertex={Z, W-minus, W-plus}; ELWK=1; FName=ZWW; 

FValue={R, " GG*GCOS"}; Vend; 
Vertex={photon, W-minus, W-plus}; ELWK=1; FName=AWW; 

FValue={R, " GE"}; Vend; 

Vertex={gluon, gluon, gluon }; QCD=1; FName=GGG; 

FValue={R, " CQCD"}; Vend; 
y 

7. vvvv 

y o 

Vertex={W-plus , W-minus, photon, photon }; ELWK=2; 

FName=WWAA ; FValue={R, " GE2"} ; Vend; 
Vertex={W-plus, W-minus, Z, photon }; ELWK=2; 

FName=WWZA; FValue={R, " GE*GG*GC0S"} ; Vend; 



y 

7. FFV (FFW : without quark mixing) 
y 



Vertex : 
Vertex^ 
Vertex^ 
Vertex : 
Vertex^ 
Vertex ; 



={positron, 

FValue={R, 
={anti-muon, 

FValue={R, 
={anti-tau, 

FValue={R, 
={nu-e-bar , 

FValue={R, 
={nu-mu-bar , 

FValue={R, 
={nu-tau-bar 



nu-e , 

GWFL"} ; 
nu-mu , 

GWFL"}; 
nu-tau, 

GWFL"}; 
electron 

GWFL"}; 
muon, 

GWFL"}; 
,tau, 



W-minus} ; 
FType="V-A' , 

W-minus} ; 
FType^'V-A' 1 

W-minus}; 
FType="V-A' 
, W-plus }; 
FType="V-A' 

W-plus }; 
FType="V-A r 

W-plus }; 



ELWK=1 ; 
; Vend; 

ELWK=1; 
; Vend ; 

ELWK=1; 
; Vend ; 

ELWK=1; 
; Vend; 

ELWK=1; 
; Vend ; 

ELWK=1 ; 



FName=WNE; 
FName=WNM; 
FName=WNT; 
FName=WEL ; 
FName=WMU; 
FName=WTA ; 
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% 

°/ Counter terms introduced by photon-Z mixing. 
% 

I SSV 

% 

Vertex={chi-3, Higgs , photon }; ELWK=3; FName=AHY; 

Vend; 

% 

•/. svv 

% 

Vertex={Higgs, Z, photon }; ELWK=3; FName=HZA; 

Vend; 

Vertex={Higgs , photon, photon }; ELWK=3; FName=HAA; 
Vend; 



% 

7. FFV (FFZ) 

% 

Vertex={nu-e-bar , nu-e, photon}; ELWK=3; FName=ANE; 

Vend; 

Vertex={nu-mu-bar , nu-mu, photon}; ELWK=3; FName=ANM; 
Vend; 

Vertex={nu-tau-bar , nu-tau, photon}; ELWK=3; FName=ANT; 
Vend; 

y======================================= 

Mend; 

y======================================= 

Lines beginning with "°/ " are comment lines. 

The file comprises of three parts: definitions of the names of the coupling 
constants, definitions of the particles and then definitions of the vertices 
and counter-terms, arranged in this order. A definition is a sequence of 
u keyword=value;" or "keyword;" . The value part may have an optional 
part, which is enclosed by brackets " [" and "]". 

The first keyword is Order, which define the names of the coupling con- 
stants. When multiple coupling constants appear in the model, a list of these 
names is specified enclosed by braces "{" and "}" . For example, 



34 



Order={QED, QCD}; 
is for a model including QED and QCD coupling constants. 

C.l Definition of particles 

The definition of a particle begins with the keyword Particle and ends with 
the keyword Pend. 

The value of the keyword Particle is the name of the particle. The 
name, beginning with an alphabet, should be unique among the particles 
defined in this file, since this name is used to identify the particle. An option 
to the name is a shorthand name of the particle, which is used to show the 
particle in the graphic output. 

In the same way, the name of an anti-particle is given as the value of the 
keyword Antiparticle. When the anti-particle coincides to the particle, 
this part should be defined as 

Antiparticle=Particle. 

The type of particle is given by the keyword PType. Its value is either 
Scalar, Vector, Majorana, Fermion or Ghost. 

The keyword Charge specifies the electric charge. The value is given as 
a signed integer or a rational number in the unit of the positron charge e. 

The dimension of color representation of the particle is specified by the 
Color keyword as an integer. 

The definition of the mass parameter of the particle comprises three items: 
the Fortran variable name of the mass, the default numerical value and a flag 
specifying massive or massless. Even if the particle is massless, a fictitious 
mass can be introduced in some part of the calculation in order to avoid an 
infrared divergence or mass singularity. We require a definition of a particle 
being massless or massive, and the Fortran variable name for the particle 
mass with a default value not only for a massive particle, but also a massless 
one. 

The keyword Massless or Massive (without a value) specifies that the 
particle is either massless or massive, respectively. If nothing is specified, the 
particle is considered to be massive. 

The Fortran name is defined by the keyword Mass and its default value 
by MValue. The value of keyword MValue is defined as a character string 
which is used in the Fortran code. 



35 



The Fortran name of the width and its default numerical value is given 
by the keywords Width and WValue, respectively. 

In the above definition, we do not assume any special name of particles; 
one can define the particle name freely. However, amplitude-generation pro- 
gram is necessary to know whether a special particle appears in a Feynman 
graph or not. For example, the CHANEL library offers a calculation of the 
amplitudes in a general covariant gauge. In the unitary gauge, x- sca l ars 
disappear form the calculation, and one must drop any Feynman graphs in- 
cluding them. In order to detect such kinds of particles, particularly in the 
amplitude generating-program, we add another keyword PCode with an in- 
teger particle code. The values of PCode are not used in the Feynman-graph 
generator in the current version. 

C.2 Definition of vertices 

A definition of a particle begins with the keyword Vertex and ends with the 
keyword Vend. 

The keyword Vertex defines the interacting particles in the list of their 
name. The direction of a particle is defined as incoming to the vertex. For 
example, 

Vertex={nu-e-bar , electron, W-plus} ; 

defines a vertex which e~ and W + are coming in and v e is going out. 

The order of the coupling constants is given by keywords defined in the 
value of the keyword Order such as: 

ELWK=1; 

The total order of coupling constant of tree 3- and 4-point vertices should be 
1 and 2, respectively. 

A program generating amplitudes requires additional information. The 
name of the coupling constant used in the Fortran code of the vertex is 
given by the FName keyword, whose data type and default value are given 
by the FValue keyword. The number of parameters in a coupling constants 
of a vertex depends on the interaction. In order to decide its number and 
how particles interact in a vertex of two fermions and a boson, the type of 
the interaction is specified by the keyword FType. The value of FType is 
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either "V", "A", "S", "P", "V-A", "V+A", "S-P", "S+P" or "NON". When this 
keyword is omitted, the fermion-boson vertex is assumed to FType="N0N". 
The format of specification by the FValue keyword is 

FValue={<i2ag>, <paraml>}; 

or 

FValue={<5ag>, <paraml>, <param2>} ; 

in accordance with the number of coupling constants. The parts <paraml> 
and <param2> are embedded into the Fortran code by the Fortran-code 
generator. The part <Aag> is either R, I or C, corresponding to the parameter 
being real, pure imaginary or complex number, respectively. These keywords 
are not used in the graph-generation program. 

One can define a counter term of renormalization as a vertex with higher 
order coupling constants. 
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D Output file 



We show an example of the output file of the graph generator for the input 
file shown in appendix |B|. 

nnnnnnnm.v.m:mm:am:mv.v.m. 

Model="all.mdl" ; 

mmmmmrammmmmmm 

Process=l ; 
External=4; 

0= initial electron; 

1= initial positron; 

2= final w-plus; 

3= final w-minus ; 
Eend; 

elwk=4;Loop=l; 
0PI=No ; Expand=Yes ; 

% 

Graph=l ; 
Sf actor=l ; 
Vertex=4; 

0={ 1 [positron] } ; 

1={ 2 [electron]}; 

2={ 3 [w-plus]}; 

3={ 4 [w-minus]}; 

4 [order={l ,0}] ={ 1 [electron] , 2[positron], 5[photon]}; 

5 [order={l ,0}] ={ 5[photon], 6[w-minus], 7[w-plus]}; 

6 [order={l ,0}] ={ 3 [w-minus] , 6 [w-plus], 8[z]}; 

7[order={l,0}]={ 4 [w-plus], 7 [w-minus], 8[z]}; 
Vend; 
Gend; 

% 

Graph=2 ; 



% 

Graph=316; 
Sf actor=l ; 
Vertex=3; 

0={ 1 [positron] }; 

1={ 2 [electron]}; 

2={ 3 [w-plus]}; 

3={ 4 [w-minus]}; 
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4 [order={l ,0}] ={ 1 [electron] , 4[w-plus], 5 [nu-e-bar] } ; 

5 [order={l ,0}] ={ 2 [positron] , 3[w-minus], 6[nu-e]}; 

6 [loop=l ; order={2 , 0>] ={ 5 [nu-e] , 6 [nu-e-bar] } ; 
Vend; 
Gend; 

% 

Pend=316; 

/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/^^ 

Process=2 ; 
External=4; 

0= initial electron; 

1= initial positron; 

2= final w-plus ; 

3= final w-minus ; 
Eend; 

elwk=2 ; Loop=0 ; 
0PI=No ; Expand=Yes ; 

% 

Graph=l ; 



Graph=28; 
Sf actor=l ; 
Vertex=2 ; 

0={ 1 [positron]}; 

1={ 2 [electron]}; 

2={ 3 [w-plus]}; 

3={ 4 [w-minus]}; 

4={ 5 [photon]}; 

5 [order={l ,0}] ={ 1 [electron] , 2 [positron], 6[z]}; 

6[order={2,0}]={ 3 [w-minus] , 4 [w-plus], 5 [photon], 6[z]}; 
Vend; 
Gend; 

% 

Pend=28; 

End=3 ; 

Lines beginning with "°/ " are ignored as comments. 
The output file is a sequence of statements. Most statements are written 
in the format 

u keyword=value" or "keyword" . 
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separated by colons " ; " . The keyword part is an identifier (name composed 
of alphabet and digits). The value part is either an identifier or a list of 
identifiers separated by "," and enclosed by braces "{" and "}". An identifier 
may have options enclosed by brackets " [" and "] " . Between two brackets, 
a statement or a list of statements is placed. 

The fist non-comment line specifies the file name of the model definition 
which is used in graph generation. 

The output file includes several processes in accordance with the input 
file. In this example, the output file includes three processes. Each of them 
begins with the line 

Process=<process number>; 

and ends with the line 

Pend=<the number of graph in the process> ; 

The description of a process comprises descriptions of external particles, 
options for graph generation copied from the input file and a sequence of 
descriptions of the generated graphs. 

The description of external particles begins with the statement 

External=< the number of external particles>; 
and ends with the statement 
Eend; 

Between them, statements for each external particle are written in the format 
<external particle number>= initial <particle name> ; 

or 

< external particle number>= final <particle name> ; 

The <external particle number> is an integer number beginning from 0. 
The description of a generated graph begins with the statement 

Graph=<grapn number> ; 

and ends with the statement 
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Gend; 

After the Graph statement, it follows the global factor of the graph, 

Sf act or=< inverse of global factor>; 

The value of Sfactor is the inverse of the symmetric factor of the graph 
with a sign corresponding to the permutation of the external fermions and 
fermion loops. 

The structure of a graph is expressed by connections between the nodes. 
They are placed between two statements: 

Vert ex=<the number of generated vertices>; 

and 

Vend; 

Corresponding to each node, information about the connections from the 
node to others is written as 

<node number>={ <internal line number> [particle name] , 

... }; 

The <node number> is sequentially numbered beginning from 0. A node 
number corresponding to one defined as <node number> expresses the same 
external particle. The <internal line number> is numbered sequentially be- 
ginning from 1. Two nodes with common <internal line number> are con- 
nected by the internal line. The <particle name> is the name of a particle 
assigned to the internal line defined as incoming to the node. When multiple 
coupling constants are defined in the model, the following option is added to 
the <node number> 

order= {order- 1, order-2, order-n}, 

which specifies the orders of coupling constants of the vertex, corresponding 
to the Order statement in the model file. Moreover, when the vertex is a 
blob or a counter term, the following option is also added: 

loop=< the number of loops in the vertex> 

The file ends with the End statement: 

End=< the number of processes> ; 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-th/9408107v2 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-th/9408107v2 



